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1 Introduction 



The large ./V expansion provides the powerful tool for the investigating of the fields theoretic 
models with a internal symmetry. Being nonperturbative in its nature this method has proved 
to be successful in revealing properties which are inaccessible in conventional perturbation theory. 
The 1/N approach relies on the expansion of the effective action in the saddle point approximation 
which determines the leading order structure in 1/N. However, beyond the lowest order the 
computations become quite untractable because the propagators of auxiliary fields depend on 
the mass in a very complicated fashion. 

To obviate these difficulties it had been suggested in Refs. Q, where the nonlinear sigma 
model had been analyzed, to consider theory directly at critical point, that should ensure the 
masslessness of all propagators. The latter circumstance simplifies considerably the computation 
of the Feynman diagrams. Solving the skeleton Dyson equation the authors of Refs. || have 
calculated the basic indices of sigma model with 1/N 2 accuracy in the arbitrary space dimensions. 
Whilst the initial application of self consistency method was to a bosonic model, the techniques 
have also been extended to examine the fermionic theories ||. The further progress in the 
calculation of the critical indices is related to the conformal bootstrap method. In the frame of 
this approach the most strong at the present moment results — 1/N 3 anomalous dimensions of 
the basic fields in the a model Q and in the Gross -Neveu model []|, ||] — had been obtained. 

If one is interested in the calculation of the critical indices related to the basic fields only, 
these methods - self - consistency equations and the conformal bootstrap method - are likely 
to be the most efficient. However, some specific features of the above methods make them 
inconvenient for the treatment of the composite operators, except for the simplest ones 0, g]. At 
the same time, the calculations of this type became actual now. First of all, it is related with 
the attempts to understand in more details the structure of the conformal fields theories (CFT) 
in D > 2 dimensions (see Refs. |9|, [H], and references therein). The basic examples of higher 
dimensional CFT being bosonic and fermionic N vector models, the 1/N expansion underlies the 
computational method used for analyzing of these models. Second, the renewing of the interest 
to the 1/N calculations is conditioned by the recent progress in the higher order computations 
of the renormalization group (RG) functions in QCD. The analytical calculations in QCD in the 



high orders (four loop f3 - function |l2j| , anomalous dimensions of twist - 2 operators [13, |14| ) 
are very complex ones because of the huge number of Feynman diagrams to be evaluated. Thus 
it is important to have methods which may ensure the independent checks of this results. The 
one of the possible approach of this type is the 1/Nf expansion Further, the necessity of 
the careful analysis of the anomalous dimensions of the composite operators of special type have 
arisen in connection with the so - called stability problem [16, 17, 15]. Recently, it was realized 
that the 1/N expansion is the more suitable approach for this purpose JlS|], 

As mentioned above, the calculations in 1/N scheme beyond the 1/N order become tractable 
only when theory is considered directly at the critical point. But in this case the corresponding 
model lost the property of the multiplicative renormalizability and cannot be worked out with 
the help of the standard RG methods. It forces one to use other approaches applicable in this 
case. The methods of self - consistency equations and conformal bootstrap, discussed above, are 
applicable for the operators of the special type only. The approach based on the operator product 
expansion used in Refs |TI| , [l0[ | is too involved to have a practical use in higher order calculations 
and limited by the models possessing of the conformal symmetry. The most convenient approach 
for the calculating of the anomalous dimensions of the composite operators in the first order of 
1/N expansion had been developed in Ref [2C[. The authors of Ref. [2C] have exploited the prop- 
erty of the scale invariance of the correlators at critical point to obtain the anomalous dimensions 
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of arbitrary composite operators. The main result of the paper [20| is the formula expressing 
1/N order anomalous dimensions via the A pole residues of the corresponding diagrams (A is 
the regularization parameter, analogous e in dimensional regularization). Certainly, this way is 
the most effective in the case of the operators of arbitrary type. 

At the straightforward generalization beyond the 1/N order (see Refs. [21, the above 
method lost many its attractive features. In the present paper we shall take advantage of the 
other approach ^] which allow to derive the simple formula for the anomalous dimensions in 
the 1/N 2 order. To avoid the inconveniences connected with the absence of the multiplicative 
renormalizability we consider extended model which is reduced to the initial one by tuning of 
some parameters. The extended model can be analyzed by the standard renormalization group 
(RG) methods. To obtain the anomalous dimensions in the initial model we relate it with the 
corresponding RG - functions of the extended one. In the following we restrict ourselves to the 
nonlinear sigma model, the generalization on the other model being straightforward. 

The paper is organized as follows: In section 2 we review the basic features of the nonlinear 
sigma model and discuss the problem to be solved. In section 3 the basic formulae for the 
anomalous dimensions of composite operators are derived. Some technical tricks useful in higher 
order calculations are collected in section 4. The results of 1/N 2 order calculations of anomalous 
dimensions of some composite operators are presented in section 5. The conclusions are given in 
section 6. 



2 Preliminaries 

The 1/N - expansion for the massless 0(N) -nonlinear sigma model in D = 2/U-dimensions is 
derived in the standard manner from the following action: 

S = -\{d4> A ) 2 - l -M~ 2 ^K^ + \^ A ) 2 + (2.1) 

where (f> is the O(N) vector field and ^-auxiliary field; A - regularization parameter, M renor- 
malization mass. The two first term in ( p.l| ) form the free part of the action, whereas two last are 
treated as a interaction. The kernel K is determined from the requirement of the cancellation of 
the simple <j> loop insertion in ip line 

• -f -o 

and reads as: 

N f d D k 1 N H 2 (l 



K(p) = -— I - " n 0/ - — - = " » D ~ 4 = b~ l p D - A , (2.2) 

2 J (2ir) k 2 {p — k) 2 2(4vr)^ H(2fi-2y 1 V ; 

where H(x) = T(fj, — x)/T(x). 

The regularized kernel is defined as = C~ 1 (A)K(p)p 2 ^~ . The function C(A) is as- 
sumed to be regular at A = and C(A) — > 1 at A — > 0. The freedom in the choice of the form 
of C(A) will be discussed in Sec.||. The propagators of (f) and ip fields have form {G^ = K^ 1 ): 

G AB (p) = 5 AB G^(p) = 5 AB p~ 2 and G^p) = 6M 2A C(A)p- 2 ^" 2 + A ) (2.3) 

in the momentum space and 

G^x) = Ax~ 2a and G^(x) = BM 2A C(A)x~ 2 ^~ A ^ (2.4) 
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in the coordinate space. Here a = // — 1, (3 = 2 and the explicit expressions for the amplitudes 
A,B can be easily deduced from the Eq. (|2.3| ). 

The renormalization of this theory carried out in the standard way. The divergences appears 
as A - pole terms in diagrams and can be excluded by the standard 1Z - operation. The power 
counting shows that the theory with action fl2.1|) is renormalizable in the 1 /N - expansion in the 
whole range 2 < D < 4 and the renormalized action has the form: 

Sr = -y(3tf - X -M~ 2 ^K^ + ^-^ 2 + \m. (2.5) 
As distinct from the more usual field theoretic models, the theory under consideration is not 



multiplicatively renormalizable. Indeed, the renormalized action (2J3) cannot been brought into 
the form of that ( |2.1| ) by a redefinition of the fields (There is not a parameter like the coupling 
constant in ( |2.lD ). Due to nonmultiplicativity, one fails to apply RG methods for the calculation 
of critical exponents and has to seek for another approaches. One of them, which is briefly 
sketched below relies on the property of the scale invariance of the theory. The latter can be 
deduced from the equivalence of the nonlinear sigma model to the (<1> 2 ) 2 theory in D = 4 — e 
dimensions at the critical point. The more formal proof not appealing to such equivalence can 
be found in Ref. p^]. Now we shall discuss the basic points of the approach mentioned above. 



(For more details see Refs.[20, 21, |i~9[.) 

Let us consider the 1-irreducible Green function T(p) that depends on the one momentum p 
only — the inverse propagator of (f) field or the vertex function with one zero momentum. We 
prefer to work with the dimensionless object: T(p) = p~ d T(p) , d being the canonical dimension 
of the function T(p). Due to the scale invariance, the form of the renormalized Green function 
is fully determined (up to some constant factor T) by its anomalous dimension 7: 

f(p) = r.(A)\ (2. 6 ) 

The differentiation of the Eq.( |2.6D with respect M yields the equation on the anomalous dimen- 
sion 

Md M T(p) = -7f (p). (2-7) 

To calculate the left side of the Eq.fl2.7l) we notice that the all dependence from M in the 
diagrams contributing to f (p) results from the propagators of ip field (see Eq.(^)), i.e. the 
every diagram acquires the factor M 2An , where n is the number of ^-lines in a diagram. Hence, 
the formula fl2.7| ) takes the form: 

7 r(p) = -2A ]T n^i, (2.8) 
where rii is the number of '(/'-lines in diagram r\. The sum runs over all diagrams Tj for the 



renormalized Green function T(p) and, of course, the limit A — > is implied in Eq. (2. 



In the first order of 1/N expansion Eq. (2.8) becomes very simple. Indeed, taking into account 
that the first order diagrams can develop the simple A poles only, Tj = (rj)i/A + (1^)0 + and 
that f = 1 + 0(1/ N) one gets 

7 = -2^n i (ri) 1 + 0(l/iV 2 ). (2.9) 

{Ti} 

Thus for the calculation of the anomalous dimensions in the lowest order one need know the 
pole parts of the corresponding diagrams only. At first sight this property is lost in the next 
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Figure 1: The example of the rearrangement of the diagrams. The black dots denote the ver- 
tex counterterm. The symbol 1Z' stands for the standard operation of the subtracting of the 
sub divergencies from a diagram. 

order. Indeed, the evaluating of the 7 with 1/N 2 accuracy requires the knowledge of the 1/N 
order corrections to the constant f, which are determined by the finite parts of the 1/N order 
diagrams. This circumstance greatly reduces the effectiveness of the Eq. fl2.S| ) for the high order 
calculations. To avoid this difficulty one may try to pick out of the sum in the Eq. ( |2.8| ) the 
diagrams reproducing the undesirable terms on the lhs of the latter. Taking in mind the analogous 
formulae in the dimensional regularization scheme we regroup the terms in the sum tfH^) in the 
way shown on the Fig |l|. It can be checked, in these particular cases, that only the TV - terms 
contribute to the anomalous dimensions, whereas the others ensure the cancellation of the extra 
terms on the lhs of the Eq. (|2.8| ). (The 1Z' is the standard operation of the subtracting of the 
sub divergencies from diagram.) Eventually, the formula for the anomalous dimensions reads 

j = -2j2n i (K'r i ) 1 + 0(l/N 3 ). (2.10) 



Here the sum runs over all diagrams (up to 1/N 2 order) generated from the action (2.1), i.e. 
the nonrenormalized diagrams. The notation is used for the residue at the simple pole 
in Laurent expansion of a function /. The formula (2.10) gives the simple algorithm for the 
calculating of the critical indeces r\ and x defined as 

70 = r7/2 and 7^ = -<q - x (2.11) 

up to 1/N 2 order: 

rj = 2^2n i (n'T H:i ) 1 , x = "2 £ n^lZ'T^^. (2.12) 
The advantages of the formula ( 2.10[ ) for the practical calculations in the comparison with 



those (2^8) is self-evident. Our next purpose is to derive the analogous formulae in the case 
of arbitrary composite operator. However, to do this we shall take the advantage of the another 
approach. 



3 The extended model 



Let us consider the model based on the action obtained from (2.1) by introducing two new 
independent couplings u and v: 

s = -\{d4> A ) 2 - \m~ 2 ^k^> + \^ 2 + v -m. (3.1) 

In the following we will refer to it as UV model. It had been firstly introduced and investigated 
in Ref. |2^] . Obviously, the initial model ( p.lQ is recovered by the special choice of parameters 
u = v = 1. As distinct from the ordinary sigma model the propagator of tp field in the UV 
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Figure 2: The effective ifi - line 



model has the more complicated form. Indeed, at v 7^ 1 the last term in the action (3.1) does 
not ensure the exact cancellation of the simple <p loop insertion in ip line and one should sum up 
all such insertions (see Fig. |2|), that yields for the propagator of ip field: 



G (p. Uj v) = I G (p) ( i + ^ — -t A + ^^t 2A + ...]. (3.2) 



Here t = C(A)(M 2 /p 2 ). The theory is evidently renormalizable and the renormalized action 
takes form: 

Sr = -y(^) 2 " \m~ 2 ^K^ + + f W, (3-3) 

where the renormalization constants depend on the coupling u,v also ( = Zi(u,v, ..)). The 
redefinition of the fields and the couplings 

$ = Z$<5>, u = uM- 2A Z u , v = vZ v , (3.4) 

where = {0,^}, ^</> = Z^ 2 , Z^ = i^Z-f 1 and Z u = Z v = Z^ 2 , brings the action ( |3~li| ) into the 
initial form, hence, the UV model is the multiplicatively renormalizable. This implies that the 
model can be analyzed by the standard RG methods. The basic RG functions are defined as 
follows: 

7$ = MOm lnZ$, fi u = M&mu = 2Au — 2wy^, f3 v = MOmv = —2^7^,. (3.5) 

The RG equation for the one particle irreducible (1PI) Green functions reads as: 

{Md M + ftjdu - n*7*)r(pi, • • • ,Pn) = 0. (3.6) 

Here we used the shorthand notations U = {u, v}, fludu = f3 u d u + Pvd v and n$7$ = n^^+n^^. 
For the set of the composite operators {Fi} mixing under renormalization the corresponding 
equation has the matrix form: 

(\Md M + Pudu ~ n* 7 *] S ik + If) r*(p;pi, ...,p„) = 0. (3.7) 

Here Ti(p;pi, . . . ,p n ) is the n-points IP I Green function with insertion of the renormalized 
operator [Fj\n. The latter is defined as [Fi\n = QikFk and Qik is chosen to ensure finiteness of 
all Green functions of operator As usual, the summation over repeated indexes is implied. 

At last the matrix of anomalous dimensions 7_f is defined as 

lF = -Md M Z F Zp\ Zf = Q ik Z~ n *' k . (3.8) 

Here Z $ n *' fe = Z, ^ Z a J l ^' k and the multiindex = {n < p > ki n i>,k} shows the number of the <j) 
and ip fields in the monomial Fk- 

We remind, that our purpose is the calculation of the critical indices of composite operators in 
the ordinary sigma model, i.e. we are interested in the limit u, v — > 1 of the UV model. However, 
in this limit the beta functions (3 U , Pv are not zero. (Indeed, (3 U = f3 v = 7^, 7^ 0.) Thus, generally 
speaking the point u = v = 1 is not the critical point of UV model and the RG functions 7$, jp 
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calculated at u = v = 1 are not coincide with the anomalous dimensions of the fields (operators) 
in the ordinary sigma model. However, it had been shown in Ref. [^] that at the particular 
choice of the renormalization scheme (the so called scheme of subtraction at fixed momenta, we 
refer to it as the scheme I) the renormalized Green functions depend on the difference u — v only. 
In this case the anomalous term 

A = -PuduTfiu - v;p,pi, . . . ,p n )\ u=v= i (3.9) 



is dropped out of the Eq. ( |3.7D and the RG functions 7$,7_f calculated in this scheme are the 
true anomalous dimensions. The Green functions calculated in any other renormalization scheme 
are related to those T\(u — v } ...) calculated in the scheme I by a finite renormalization: 

Fi(u, v, ...) = Z* k (u, v) • T{(u - v, ...). (3.10) 

Here u = Z*(u,v) ■ u , v = Z*(u,v) ■ v and Z* k (u,v), Z*(u,v), Z*(u,v) are the constants of 
the finite renormalization. Taking into account that the bare couplings uq,vq are independent 



from the renormalization scheme, we obtain from the Eq. (3.4) that u/v = u/v and, hence 



Z* = Z*. Then, from the Eq. (|31C| ), one immediately obtains that Green functions calculated 



in the different regularization schemes at u = v = 1 differ one from another only by a constant 
multiplier. This proves that the critical exponents calculated in the model based on action (|2.1| ) 
are independent from the regularization scheme being used and, in particularly, from the choice 
of the function C(A). 

Further, inserting Tj in the form ( |3.10| ) into Eq. ( |3.7| ) one obtains the following equation for 
the Green functions of the ordinary sigma model: 

Md M ~ rw] $ ik + (7f + IT)) r fe (p;Pl, • • • ,Pn) = o. (3.11) 

where 7^ = {f3 u d u + Pvd v )Z*(u, v)Z*~ 1 (u, v)\ u=v= \. Thus the anomalous dimensions are given 
by the eigenvalues of the matrix 7^ + 7^, but not of the 7^ alone. Of course, the 7^ is zero 
in the scheme I, but the latter is not very suitable for the practical calculations. The most 
convenient scheme for the calculations of RG functions is the MS scheme, but then one has to 
know -fp. Fortunately, the matrix 7^ is appeared to be of 0(l/iV 3 ) in the MS scheme, i.e. for 
the determination of the critical exponents up to 1/N 2 order it is sufficient to know 7^ only. 
Henceforth we restrict ourselves to MS scheme only. The RG functions 7, then, are related to 
the simple A poles of the corresponding renormalization constants: 



^ = 2d u zW ^=2d u Z^ 1 f=-2d u Q^ +^7*. (3.12) 

T U=V=1 v U=V=\ U=V=\ 

where Zfav) = 1 + En=i zl n \u,v)/A n , Q lk (u,v) = 5 ik + En=i (u, v) / A n . The mixing 
matrix is calculated directly from the diagrams as the divergent part of the functional 



AD 



00 -t 

K.n'T Fi (x; $) = - £ — • Q$F k {x), (3.13) 

n=l 

F Fi (x;$) = X!~T [ dxi..dx n r Fitn (xixi...x n )$(x{)...®(x n ). 



Here T Fi (x;xi, . . . ,x n ) is the nonrenormalized (i.e. one arising at averaging with action ( |3.1[) ) 
n-points 1PI Green function with operator Fi insertion and the operation /C selects the singular 
(A poles) part of the diagrams. 
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Because there are not derivatives with respect v in the formulae ( p. 12 ), we can put v = 1 



from the very beginning. In this case the propagator of if) field (|3.2|) is reduced up to factor 
1/u to that in the ordinary sigma model: G^(p,u,v = 1) = G^(p)/u. Then taking into account 
that the constants Zi,Qik are determined by the divergent (pole) parts of the corresponding 
Green functions and that operation d u gives simply the factor — (nwj is the number of ip 
lines) for each diagram, one concludes that all calculations needed for determining of anomalous 
dimensions ( 3.12j ) can be done in the framework of the standard a model. 



Let us prove now that matrix 7^ vanishes in the 1/N 2 order. For this it is sufficient to show 
that the anomalous term A (see Eq. (|3.9| )) is zero in this order. First of all, note, 7^ being of 
0(1/N) order, we need to prove that the remaining part is of 0{1/N 2 ) order. Let us consider 
the arbitrary first order diagram of the ordinary a model with n internal ifi lines. Regarding 
the regulator A on each line as the independent variable and picking out the multiplier C(A) 
from the i/j propagator we write the answer for the diagram as C n (A)G(A, . . . , A). Then, taking 
into account the form of the ip propagator Eq. ([T^) one obtains the following expression for this 
diagram in the extended model: 

1 00 1 

— Y (C(A)) mi+ - +m "G(miA,...,m n A)(^^) mi+ - +m "- n . (3.14) 
u n . u 

mi,...,m„=i 

As usual, the renormalized diagram is obtained from the nonrenormalized one by the 1Z operation. 



In the Ref. [22] it was shown that the 1Z operation based on the subtraction scheme I gives for 



each term in the sum ( |3.14j ) the same (up to factor {{v — l)/n) ") answer. (Then, the sum ( |3.14 ) 



becomes a trivial and results in the factor (1 + u — v) n .) The above property holds no longer in 
the MS scheme. The action of the 1Z operation on different terms in the sum (gig ) leads to the 



essentially different answer even in the first order of 1/N expansion. However, in the first order 



of the 1/N expansion the whole sum (3.14) after renormalization depends on the difference u — v 
only too. 

In the case of a convergent diagram G this statement is trivial. Let us consider the superficially 
divergent diagram. In the 1/N order such diagram having no subdivergencies, we represent 
function G(. . .) as: 

fit a A \ F(miA, . . . ,m n A) 

G(miA, ...,m n A) = ■ ■ — -, (3.15) 

(mi + . . . + m n )A 

where the function F is regular in the neighborhood of zero in all its arguments. Since in the MS 
scheme the R operation is reduced to the removing of the A poles one has for the renormalized 
diagram: 

^ mi ,.i„=A m 1 + ... + m n J { u 

Taking into account that the terms like m\/ (jn\ + . . . + m n ) due to the evident symmetry of the 
rest of summand in rrij can be replaced by 1/n, one obtains for the sum (|3.16|) 



nC'(0)F(0)+F{(0) + ... + Fl l (0) 

n{l + u-v) n 1 ' ' 

The case when a superficially convergent diagram has a divergent subgraph can be treated along 
these lines as well. 

Since the anomalous term A is proportional to (d u + d v )T{pi, u, v) \ u =v=i we conclude that A, 
and hence, 7* vanish up to 1/N 2 order. (Note also, that we proved the more strong statement 
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Figure 3: The 1/N - diagrams for the Green function with the insertion of the operator <\>a4>b- 



A(u,v) = at u = v , than it was needed: A(u,v) = at u = v = 1.) Unfortunately, the 
arguments used above do not work in the next orders of 1/N expansion, since the diagrams 
which give the nonzero contributions to the anomalous term A can be easy found in the higher 
orders. Though we cannot to reject the possibility of the cancellation of the contributions from 
the different diagrams, it seems to be very unlikely. 

Thus we have shown that up to 1/N 2 order the anomalous dimensions of composite oper- 
ators are simply related to the A poles of the corresponding diagrams. The simplicity of the 
formulae (3.12) in the combination with the uniqueness of the triple vertex (at A = 0) and 
masslessness of propagators makes the 1/N 2 order calculations a rather straightforward task. 
The only thing which causes some problem is a large number of diagrams to be evaluated. In 
the next section we show how to overcome this difficulty. 



4 Conformal structures in perturbation theory. 

Everyone being a familiar with the 1/N expansion knows that the number of diagrams to be 
evaluated increases drastically with a order of expansion. To receive the impression about the 
rate of the growth it is sufficient to compare the number of the first and the second order diagram 
for any operator. Indeed, one has only one 1/N order diagram (Fig. ||) and eight ones in 1/N 2 
order (Fig. |^) for the 4> A (f) B operator. The analogous comparison for the tp 2 operator gives 3 
and 69 diagrams, respectively. However, it can be easy realized that the bulk of the second 
order diagrams are the nothing else as the first order ones with the insertion of the vertex or the 
propagator corrections. It is intuitively clear that taking into account of these corrections must 
lead to the dressing of the bare propagators and vertices. In what follows we shall show that 
this really occurs. For further discussion to be more transparent we will use (f> A <j) B operator as 
illustrative example. The second order diagrams for the two-point Green function with operator 
cj) A (j) B insertion are drawn on Fig. ||. Two of them (c,d) arise due to the vertex corrections to 
the first order diagram (Fig. |3[); one (e) - due to the self-energy (SE) insertion in <j> line and 
three (f,g,h) - due to the SE insertions in ifi line. 

We are interested in the contributions from the diagrams of this type (c - g on Fig. ||) to the 
matrix d u Q lk \ u= \. The more precisely we need the simple A pole term in the sum 



Here sum runs over the 1/N 2 order nonrenormalized diagrams {T{} which arise from the 1/N 
order - base - diagrams after taken into account the vertex and the propagator corrections (the 
diagrams c-h on the Fig. ||). For later convenience we rewrite the sum (4.1) in the following form: 



E 




(4.2) 
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(a) 



(b) 



(d) 






(e) (f) (g) (h) 

Figure 4: The 1/N 2 - diagrams for the Green function with insertion of operator 



Here {f j} is the set of the counterterm diagrams corresponding to the ones from the set {1^}, i.e 
IZ'Ti = Ti — Tf, fii is the number of tp lines in the diagram Fj, while hi = rii — hi is the number 
of ip lines in the divergent subgraph. It is appeared that the last sum in the Eq. ( |4,2| ) has no 
simple A pole, provided that the relation 2n^ = = ny between the number of the internal 
lines and the triple vertices holds for the 1/N order base diagram. Indeed, one can see that the 
divergent subgraphs in the diagrams with the SE insertions in tp line (Fig. ^| f-h) are the same 
as those in the diagrams with the vertex corrections and the SE insertion in (f> line (Fig. || c-e). 
We combine the diagrams into three different groups — (c,f), (d,g), (e,h) - in accordance with 
the type of the divergent subgraph. Obviously, the diagrams appearing after reducing of the 
divergent subgraph in the diagrams from the same group, i.e. those from the set {Ti}, differ one 
from another only by the simple (j) loop insertion in ip line. For example, for S/A being the pole 
part of the divergent subgraph, the contribution of the (e,h) pair into the sum in question reads: 





+ 



(The factor 1/2 here is due to the symmetry coefficient of diagram (h) 
expression for the above sum can be represented as 



The analytical 



(4.3) 



F being the regular function of A. (To deriv e t}4.3|) it is sufficient to remember the defini- 
tions (2.2),(2.3), (|2.4| )). It is seen the difference (|4.3|) does not contain the simple A pole. Obvi- 
ously, the same arguments work and in a general case, the condition 2n^ = = ny ensuring the 
cancellation of the simple pole contributions among each group separately. Henceforth we shall 
only deal with the diagrams satisfying this condition. Moreover, we confine our consideration by 
those diagrams for which the only sub divergencies are due to the vertex and SE subgraphs. In 
spite of that these conditions restrict the class of the admissible operators it appears to be wide 
enough to justify the necessity of the present discussion. 

We start with the analyzing of the diagrams with the SE insertions, say, in (j) lines (Fig. ||] e). 
For later convenience we will regard the line indices as the variables and distinguish the regulators 
A on the ip - lines being internal (Ai) and external (A2) with respect to the divergent subgraph. 
Also, we introduce the regular function F related to the base diagram (Fig. ||) with the line 



10 



(a) (b) 
Figure 5: The 1/N order vertex diagrams 



indeces differing from their canonical values au = 1 and = fi — 2 by Aj as follows 

q = F{a<i> + Ai, a + A 2 , aty + A 3 ) 
A1 + A2 + A3 

Then the answer for the diagram with the SE insertion in <j) line can be written as 

5(AQ F{ a<j) + Ai,q q v , + A 2 ) 
Ai A 1 + A 2 ' l4 - Dj 

where the factor 5(Ai)/Ai result from the integration over loop momenta in the subgraph. 
For the corresponding counterterm diagram one has 5(0)/Ai • F(a^, a^a^ + A2VA2. Their 
contribution to the sum ( [4 .2] ) reads (n = 2 and n = 1 in this case) 



2 . ,g ( Al ) . -F(o^ + Ai,Q0,q^, + A 2 ) _ 5(0) _ F{a^a^a^ + A 2 ) 
Ai Ai + A 2 ' Ai A 2 

Evaluating of the above difference at the "symmetric" point Ai = A 2 = A yields for the simple 
pole residue 

S{^)F[{a^a^a^) + S f (0)F(a ( f > , a^, a^). (4.7) 

Noticing that the constant —5(0) is nothing else as the anomalous dimension of (f> field: 7^ = 
— 5(0) + 0{1/N 2 ) we relate the above expression to the base diagram with one bare 4> propagator 
being replaced by the dressed one: 1/p 2 -» A/p 2 ^ (A = l + 5'(0) and 7/ = 2^). Namely, the 
expression (4.7) can be obtained as the 1/N 2 - order term in the expansion of the following object 



A ■ F(a ( j > - 70,a0,a^). (4.8) 

It is not more hard to see that taking into account the SE insertions in ip line also leads to 
the dressing of this line, i.e. the contributions of these diagrams to the sum ( |4.2[ ) can be received 
from the expansion of the object: 

B ■ F(a 4> ,a 4n a^ - 7^,). (4.9) 

The following necessary step is the inclusion into consideration of the vertex corrections. 
Indeed, having restricted by the propagators corrections only one loses the more attractive feature 
of 1 /N - expansion — the uniqueness of the triple vertex. The analysis of the vertex corrections 
is a bit more complicated than those of the SE insertions. Indeed, the diagram with insertion 
in a line has, in fact, the same topology as the initial one, and hence, (that is crucial) they both 
can be described by one and the same function. To achieve this in the case in question, where 
the topology of the diagrams is essentially different, we proceed in the following way. First of 
all, let us consider the vertex function itself. Up to terms of 0(1/N 2 ) order it reads 

/ \ -1 1 7i(A;p,g) 72(A;p,q) , 1 , 

7fl(P,9) = 1 + ^ + 2A P = l + ll ^ p ' q > + 2 l2 ^ p,q '- ( °) 
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Here 7i(A)/A, 72(A)/2A are the contributions of the diagrams shown on the Fig. [|a,b, respec- 
tively. (In the further we shall not show explicitly the momentum dependence in the functions 
7i,7ii. Note only, that at A = the functions 7, are independent from p,q.) Due to conformal 
invariance of the theory, the form of the full vertex 7^ is fixed up to the constant factor. In 
coordinate space it reads 

lR (x, y; z) = Z{z - x)- 2a (z - y)- 2a (x - y)~ 2 ^ (4.11) 

where a = p, - 1 - 7 ^/2, (3 = 2 - 7^ + 7^/2. 

Further, taking in mind the pole structure of the diagrams (c,d) on Fig. || we write for the 
latter 



F[ 7m (Ai);A 2 ] 
mAi(mAi + A 2 ) 



(4.12) 



where m = 1,2, respectively. The function F is regular with respect the variables A 1 ( 2 ) and is 
defined as follows : 

F[ 7m (Ai);A 2 ] = (mAi + A 2 ) J -^- lm (A h p,p + ?)(G/7i)(p, A 2 ). (4.13) 

Here the function (G/^\){p, ....A 2 ) is related to the diagram with the vertex subgraph removed. 
(For the diagram on Fig. ||c,d it is given by the product of three propagators). The evaluation 
of the contribution of diagrams (c,d) on Fig. together with the corresponding counterterm 



diagrams yields for the simple pole residue in the sum (4.2) 



+ A 2 | ^-^(p,p + q )(G/ 7l )(p,...,A 2 ), (4.15) 



d A (F[7i(A);0] + 1% 2 (A);0]) | a=q . (4.14) 

Looking on the Eq. (|4.13| ) one can see that the renormalized vertex 7^ should appear in the 
natural way if it were be possible to carry out differentiating with respect A under the sign of 
integral. However, at A 2 = the integral evolves the pole in Ai, therefore this operation is not 
allowed. To obviate this difficulty, one can, taking into account the regularity of the function F 
in A]f2)> carry out differentiation at the nonzero second argument A 2 and then take the limit 
A 2 — »• 0. At first step we obtain for ( |4.14[ ) 

%i(0)+72(0),A 2 ) f d 2 »p 
A 2 + A2 J (2^ 

where 7^ = 7i + 7 2 /2 (see Eq. ( |4.1C| )). To proceed further we represent 7 jjP as 7^ = 
dtlR{t/N)\t=o and, again, change the order of the differentiation and integration. At last, taking 
into account that scaling dimension x ( 2.12 ) of the function tr is equal to —2(71(0) + 72(0)) at 
first 1/N order, so that the integral in Q4.1E ) evolves the pole in (A 2 — tx/%), we find eventually 
for ( [4.14| ) the following representation 

dt (-f / ^7«P l P,P + «)(G/7i)(p 1 ..,0)j [ =0 - (4.16) 

In other words the contribution of the diagrams in question is obtained from the 1/N expansion 
of the diagram with the dressed vertex. 

Basing on these results (Eqs. ([!.£]), (4. 9), ( |4. 16] )) we can formulate the simple algorithm for the 



calculating of the contribution of the diagrams with the propagator and vertex corrections to the 
matrix of the anomalous dimensions. Let G to be a 1/N order diagram and the diagrams {Gi} 
are the 1/A^ 2 order ones from the class under consideration. Then to obtain the contribution of 
the diagrams {Gi} to the anomalous dimension one need: 
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Figure 6: The confer mal triangle. 



1 . replace in the base diagram G the bare propagators by dressed ones and the pointlike triple 
vertex by the conformal triangle (see Fig. ||): 

A ~ B 

G ^ X ) = ^2A^> G *( x ) = ( 4 - 17 ) 

lR (x, y; z) = Z(z - x)- 2a {z - y y 2a {x - y)~ 2li . (4.18) 

Here a = (jl — 1 — 7^/2, (3 = 2 — 7^ + 7^/2 and A^, A^ are the full scaling dimensions of 
the fields <j>, ip, respectively. A^ = \i — 1 + 7^, A^ = 2 + 7^. 

2. introduce a regulator A in anyone line (It is needed because after the above changes one 
obtains a superficially divergent diagram). 

3. evaluate the A pole residue in the resulting diagram and pick out the 1/N 2 order term - 
(j( 2 ) . Then the required contribution is given by -2G( 2 ). 

The explicit expressions for the amplitudes A , B , Z are given in Appendix A. 

Following these prescriptions one need calculate only one — " conformal" diagram (see Fig. [?]) - 
instead of the six 1/N 2 order diagrams shown on Fig. [| (c-h). It should be noted also, that the 
resulting diagram, unlike the original ones, has superficial divergence only, that, in its turn, gives 
the definite advantages at the calculations. We end this section with remark, that the above 
procedure being useful even if the condition 2n^ = = ny is not fulfilled. Indeed, the only 



changes which will occur are that the last sum in the Eq. (4.2) will not be zero any longer and 



will contribute to the anomalous dimension on the equal foot with the first ones. 



— * —* — * — * 

5 The anomalous dimensions of the (®$) s and ($ <g> {<g>d) n Q) op- 
erators in 1/N 2 order 

In this section we apply the technique developed above to calculate the anomalous dimensions 
of the operators 

Ff = (®c£>) s ee Sym0 Al ...<f) As - traces, (5.19) 
F 2 n = ($ ® (<g>a) n $) ee Symc^^ . . . d in 4> B - traces (5.20) 

in the \/N 2 order. Here, the symbol Sym implies the symmetrization over the internal indeces 
as well as the spatial ones. 

The operators F\ are multiplicatively renormalized, whereas F2 does not. However, the 
operators admixing to F2 being the total derivatives, they are irrelevant for our purposes. The 
diagrams to be evaluated for the determination of the anomalous dimensions of operators F\, F2 
are those shown on Fig. || and Fig. ||, the latter being obviously relevant only for the Fx operators. 
Because the diagrams c — h on Fig. |] can be replaced by one "conformal" diagram (Fig |7|) in 
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A«i + A 



Figure 7: The " conformal" diagram. The black circles denote the conformal triangle. 

accordance with the rules formulated in the previous section, one has to calculate 5 diagrams for 
the F\ operators and 3 ones for the Fi operators. 

We do not dwell upon the technique of the calculating of the diagrams, since it is discussed in 
very details in the papers Jl], ||, ||, and give the final results for the anomalous dimensions 

only. (For completeness, we give the explicit expressions for the pole terms of each diagram in 
Appendix [l| and the some useful formulae are collected in Appendix 

The anomalous dimensions 7( s ) and 7( n ) of the operators F( and F% it is convenient to 
represent as 

7(s) = u (s) + sjcf, 7(n) = u (n) + 270, (5-21) 

where 7^ is the anomalous dimension of the field <f> known up to 1 /N 3 order [|| . Then representing 
u in the following form u = Ylh=i u ^ /N k we have for the two first terms: 

■4) = ~w^^ (5 - 22) 

«W = -^ 4 ( /-V(^-2)^ 2(S " 2W/(1) "^ (/i)) 

+ [n(2fi _ 3) + 2(m - 1) (2fi 2 - 3/z + 2)12] } , (5.23) 



U L = -^17 1 T^w - o I — \\— " o ' ~' R-— — ~~~~T7" — —S(n,n) 



>) 



M(m-1) r(2 M 2 -3// + 2)_ 2(/i-l)(2 / u-l) 



(/i - 1 + n)(// - 2 + n) I /i-2 ' fj, - 2 

^_ 1)B( „, M)+ ^(_l_ + *^±i; 

K/x-l) ^ 1 1 



2(/i - 1 + n)(/i - 2 + n) \ (// - 1 + n) (/x - 2 + n) 



)}, (5-25) 



where 



m = 4(2 - M )r(2/i - 2)/r 2 (/i - i)r(2 - //)r( M + 1), (5.26) 

R = + i>(n - 1) - V(2 - n) - ip(2fi - 2), (5.27) 
S(n,n)=ip(fi-l + n)-tp(n-l), (5.28) 



1 /-i 



R(n,n)= / dad/3a^- d /3^(l-a-/3) n . (5.29) 
jo jo 

The critical exponents of the operators F\, F2 are known with the high accuracy in the 4 — e 
as well as in the 2 + e expansions. Expanding our result for jr s \ in the e series near d = 2 we 
have obtained the full agreement with the four - loop result of Wegner [25]. The same exponents 
are known up to e 3 order in the 4 — e expansion [^(| . In this case we find out the discrepancy of 
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(a) (b) 
Figure 8: The 1/N 2 order diagrams with insertion of operator F\. 



our result with those of Wallace and Zia in term ~ e 3 s(s — 1)/N 2 (in our answer this term 
enters with the coefficient (—10), against the coefficient (—8) in the paper [pq] ). 

Further, using the results of the paper [27|, where the four - loop anomalous dimensions 
of 4>(§§d) n 4> operators had been calculated in the scalar cf) A model in the 4 — e expansion, we 
reconstructed the answer for the 0{N) vector model by restoring the O(N) coefficients for all 
diagrams. Comparing the expressions for the critical exponents of F2 operators derived in the 
both (1/N and 4 — e) approaches we found that they are in the full agreement up to e 4 /iV 2 terms. 



6 Conclusion 

In the present paper we have developed the simple method for the calculating of the anomalous 
dimensions (up to 1/N 2 order) in the 1/N expansion. In the comparison with the other well 
known methods of the conformal bootstrap and the selfconsistency equations this one has the 
advantage of being not restricted to the special class of the operators. In this approach the 
formulae relating the anomalous dimensions with the divergencies of the corresponding Feynman 
diagrams are very simple and in the one to one correspondence with their counterparts in the 
dimensional regularization scheme. These features together with the uniqueness of the triple 
vertex of the nonlinear sigma model make the 1/N 2 order calculations of the rather straight- 
forward task. In some cases (see for example the sec. ||) the calculations are not much more 
complicated than the two - loop ones in the O(N) model in the 4 — e expansion. However, in a 
general situation, beyond 1/N order, the number of the diagrams to be evaluated is usually large 
enough. We have shown that this problem is partially removed after going on to the diagrams 
with the dressed propagators and vertices. The number of the diagrams is considerably reduced 
by this trick and only slightly exceeds the number of those in the conformal bootstrap method, 
which is considered to be the most economic in this sense. 

We conclude with remark, that though so far we consider the sigma model only, it is evident 
that the same technique can be apply to the many other models admitting 1/N expansion - 
CP(N), Gross - Neveu and Thirring models, etc. 
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Figure 9: The 1/N order diagrams for the Green functions , and Y^^. 



Appendices 



A The renormalized Green functions in 1/N order 



In this appendix we give the explicit expressions for the amplitudes A, B, Z introduced into 
consideration in the Sec. |j. Obviously, these quantities depend on the our choice of the function 
C(A) (see Eqs. ( |2.3| ),(^4|)), However, it is easy to understand, that in the first 1/N order this 
dependence is almost trivial and can be taking into account by the redefinition of the parameter 
M. We specify C(A) by the condition for the propagator of ip field in coordinate space (Eq.( |2.4D ) 
to have the amplitude independent from A, except from the factor M 2A . Then one has for the 
amplitude A, B of the bare propagators in the coordinate space: 



A 



HQ) 
4ttM ' 



B 



-M 



2A 32ff(2-/i) 
H(2)H 2 (1) ' 



Evaluating the 1/N order diagrams for IP I Green functions T 9 
we obtain for the amplitudes A B Z: 



(A.l) 

shown on Fig. | 



A 



B 



AM~ 2 ^ ■ 
BM~ 2 ^ 



1 



7/1 fJ. 2 + fJL - 1 



2N n(n - 1) 



1 + 



7/1 / 2fi 2 - 3/i + 2 



r2yX g 2 (i)%-2)rW 

2 7T 2 ^ 



■ R 
1 + 



4/i 5 - 19/i 4 + 25/i 3 -\i 2 - 14/x + 4 X 
M^-l)(^-2) 2 / 
7/1 ( / u-3)(6/j 2 -9// + 2)" 



N 



Oi - 2)= 



(A.2) 
,(A.3) 
(A.4) 



Here 7/i, i? are given by Eq. ( |5.26|) ,( 5l27| ) and 7^, 7^ are the known anomalous dimensions of the 
fields 4>, tp 



14> 
7V> 



7^ + 0(l/iV 2 ), 
2r/i (M-1)(2M-1) 



+ 0(l/iV 2 ). 



iV 2-/7 
The index x is defined as x = ~~ (7i/> + 27^). 



(A.5) 



B Values for the 1/N 2 order integrals 

In this appendix we give the explicit expressions for the contributions to the anomalous dimen- 
sions 7( s ),7( n ) of the corresponding Feynman integrals. We will use the labels A, B,C, D, E for 
the diagrams on the Fig. |] a,b, Fig. fj] and Fig. |8] a,b, resp., as well as for the values of the 
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(s) 

latters. In the case of the F± operators we have obtained for the simple A pole residue of the 
corresponding integrals in the units rj 2 /N 2 : 



A 



B 



C 



D 



E 



s{s- 1) /i 2 (/i 2 -5/i + 5) 
2 8Gu-l)(/i-2)3' 

8(3-1)^-1) 



8(^-2)3 



8(8-1) 



(2/i -2)2 
s(s - l)(s - 2) 



/i(/i 2 - 5/x + 5) 
(/x-l)(/i-2) 

^ 2 



+ (2/x 2 - 3/x + 2)i? 



8(/i-l)(/i-2) 2 ' 



a(a - l)(g - 2) 3/i 2 (/i - i; 
6 2(/i-2) 2 



C. 



(n) 

For diagrams with the insertion of the F% operators we derived 



A 



B 



C 



m 2 (m-i) 



8(/i-l + n) 2 (/i-2 + n) 2 
/i 2 (/i-l) 2 



H — l + n /i — 2 + n 



R(n;n), 



8(/i- l + n)(fi-2 + n 

li(n - 1) r (2/i 2 - 3/i + 2) 



(B.l) 
(B.2) 
(B.3) 
(B.4) 
(B.5) 

(B.6) 
(B.7) 



2(/i- l + n)Qu-2 + rc 
2(/i-l)(2 M -l) 



/i-2 



S(n,/i) + 



/x-2 
1 / 



R 



Hifj 2 - 5// + 5) 
/i-2 U-2 + n ' - l)(/i - 2) 



+ 



}• 



(B.8) 



Each of these quantities enters into expression for the anomalous dimension with the coefficient 
(— 2n^), if we accept ru, = 1 for diagram C. 

C Chain rule and uniqueness relation 

Here we collect some rules (see Refs. [||, |J ^]) which form the basis for the calculations of 
diagrams in the massless theories. As usual, we will denote the propagator l/\x — y\ 2a by the 
line with the index a. 

The Fourier-transform of |x| _2a reads 



e i P x n 2 2a 'H(a) 



2a 



P 



,2a' 



(C.9) 



We remind that \i = D/2, H(z) = T(/i — z)/T(z) and z' = (jl — z. 

Further, we consider the chain of propagators, which is illustrated on Fig. 10. The point where 
the ends of propagators are joined is the vertex of integration. The result of this integration is 
a new line with index a + (3 — \i and amplitude proportional to 7r^i7(a, (3, 2/i — a — f3), where 
H(ai,a 2 , a 3 ) = rjf=i 
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a 



3_ 



a + (3 — H 



7 



tt^H(cx, (3, 7) 0' W 



7 

a + /3 + 7 = 2// 
Figure 10: Chain rule and uniqueness relation 



The other useful relation is the so-called uniqueness relation. In the case if the sum of the 
exponents of the lines meeting at the vertex is equal to the dimension of space-time it can be 
related to a triangle graph (see Fig. 10). 

The corresponding generalization of these rules needed for the evaluation of the the diagrams 
with the insertion of the tensor operators, are given in Refs. [24]. We give here only that which is 
relevant for the treatment of the traceless and symmetric operators, like operator F2 in our case. 
It is convenient to swamp the the spatial indices v\, v n of a operator by contracting with a con- 
stant vector u l , such that u 2 = 0, and consider scalar operator T n {x,u) = T^ 1 '" Un (x)u^ 1 . . - u Vn . 
In the course of calculations one should often to integrate the chains of the propagators like 
(ux) k /\x\ 2a . Then, the corresponding rule for the integration of the chains reads 



(x, u) n (z — x, u)' r 



x 2a (x 



)23 



7T 



, H n (a)H m ((3) {u,z) m + n 
H n+m (a + p-fi) z 2 («+^) 



(CIO) 



where H n (z) = T(/j, + n — z) /T(z). 
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